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Abstract. We observed a field-induced staggered magnetization in the 2D frustrated 
dimcr sing let spin system SrCu2(B03)2 by "B NMR, from which the magnitudes 
of the intradimer Dzyaloshinsky-Moriya interaction and the staggered (7-tensor were 
determined. These anisotropic interactions cause singlet-triplet mixing and eliminate 
a quantum phase transition at the expected critical field He for gap closing. They 
provide a quantitative account for some puzzling phenomena such as the onset of a 
uniform magnetization below He and the persistence of the excitation gap above He- 
The gap was accurately determined from the activation energy of the nuclear relaxation 
rate. 

Spin systems with singlet ground states exhibit a variety of quantum phase 
transitions in magnetic field p. A generic example is the Bose- Einstein condensation of 
triplets when the field exceeds the critical value at which the excitation energy vanishes. 
This results in an antiferromagnetic order with the staggered moment perpendicular to 
the field, as has been observed, e.g., in TlCuCla Another possibility is the formation 
of a superlattice of localized triplets due to repulsive interactions, which translates into 
magnetization plateaus at fractional values of the saturated magnetization. The best 
known example is SrCu2(B03)2 with its two-dimensional network of orthogonal dimers 
of Cu^+ ions (spin 1/2). This material shows an excitation gap Ao=35 K and plateaus 
at 1/8, 1/4, and 1/3 of the saturated magnetization [SlIllEI- A magnetic superlattice 
at the 1/8-plateau has actually been observed by NMR experiments jB]. 
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Figure 1. (color online) The magnetic layer of SrCu2 (603)2 viewed along (a) the 
c- (z-) and (b) the [110]- (y-) directions. Numbers are the site indices for both Cu 
and B. Symbols for D and D' indicate the direction of d in the Dzyaloshinsky-Moriya 
interaction d-{siXSj) with the bond direction i — > j shown by arrows. 



These basic properties can be explained by a Heisenberg model on the frustrated 
Shastry-Sutherland lattice, 

Ho = J^s,- Sj + J' ^ Si - Sj, (1) 

n.n. n.n.n. 

where J and J' are the intra- and interdimer exchange interactions as shown in figure 
jHl n\- However, several aspects of low-temperature and high-field properties remain 
mysterious: i) a finite magnetization appears well below the expected critical field for 
the gap closing = Aq/ g^fis, where is the 5^- value along the field direction |3E|; h) 
a gap seems to persist above He UHl IH] ; iii) there appears no phase transition down 
to T=0 at fields below the 1/8-plateau, which means that no Bose condensation occurs, 
while one has been observed above it JT]; iv) the magnetization shows a discontinuous 
jump at the lower boundary of the 1/8-plateau [01 Ej- The properties i) and ii) suggest 
that triplet states are mixed into the ground state by some anisotropic interactions. The 
interdimer Dzyaloshinsky-Moriya (DM) interaction, 

(2) 

n.n.n. 

where | D'- \= D' and the sign alternates as shown in figure^, was required to explain 
the splitting of triplet energy levels ^21 IH] • However, it does not have matrix elements 
between singlet and one-triplet states. Although the intradimer DM interaction is a 
candidate, no estimate for its magnitude has been made so far. 

In this Letter, we report observation of a field-induced staggered magnetization 
by ^^B NMR experiments, which is also caused by singlet-triplet mixing. Quantitative 
estimates of the intradimer DM interaction and the staggered (7-tensor are obtained. 
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Furthermore, these anisotropic interactions provide quantitative accounts for the above 
points i) and ii). We also report the behavior of the gap obtained from the nuclear 
relaxation rate upon entering into the 1/8-plateau. 

A single crystal of SrCu2 (603)2 was grown by the traveling-solvent-floating-zone 
method using LiB02 solvent The NMR measurements below 18 T were performed 
at ISSP, University of Tokyo, while data at higher fields were obtained using a 20MW 
resistive magnet at the Grenoble High Magnetic Field Laboratory. 

We first describe the anisotropic interactions compatible with the crystal structure 
(space group IA2m [14j). A crucial feature is the non-coplanar buckling of the magnetic 
CuBOs layers as depicted in figure^. This allows additional anisotropic interactions, 

4 

'Hi= - hbH ■ (^ gi ■ Si) + 

i=l 

D |E(^i^2 - s^s^,) - J2{slsl - s^sl)^ , (3) 

where A (B) denotes dimers along the x- (y) direction. The second term, the intradimer 
DM interaction, is allowed since the buckling breaks inversion symmetry at the center 
of dimer bonds. The first, Zeeman term involves anisotropic ^-tensors, which for site 1 
is given by 

^ gx -gs ^ 
gi = gy 

\ -gs gz J 

while g2, gs, and g4 are obtained from gi by symmetry operations of the crystal 
(figure Q). The sign of the x2;-component is opposite for g2, i.e., gs represents the 
staggered component. The diagonal components were determined as gx=gy=2.05 and 
gz=2.28 from ESR measurements jH]. We estimate (78=0.023 by assuming that the 
principal axis of gi coincides with that of the electric field gradient at Cul nuclei (5.6° 
from the c-axis |15J), which is nearly perpendicular to the plane containing the four 
oxygen atoms surrounding one Cul. 

In magnetic fields, the anisotropic Hamiltonian Tii results in a non-coUinear 
magnetization consisting of four sublattices mj=gj ■ {si)fiB {i=l to 4). We define 
the uniform and staggered moments on A and B dimers as = (mi -|- m2)/2, 
mf = (mi — m2)/2 and = (m^ + m4)/2, mf = (m^ — 1714) /2. The crystal 
symmetry tells us that an external field i/cxt || z (iJcxt || x) induces staggered moments 

II X and mf \\ y {mf \\ z). Such non-collinear moments, unobservable from bulk 
magnetization measurements, can be accurately detected by ^^B NMR. The magnetic 
hyperfine field acting on a nuclear spin at site i is expressed as Hf^=Y^^j=-^^ij ■ mj, 
where Aj^ is the hyperfine coupling tensor between the i-th nuclear spin and m^. The 
projection of H^^ along iJext, denoted as Hi, is obtained from the shift of the NMR 
frequency. 
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Figure 2. Angular variation of Hi at different temperatures in the external field of 
6.948 T. 



We first discuss the results for a small field (6.948 T) rotated in the (110)- (xz-) 
plane. The NMR spectra consist of several sets of quadrupole-split lines. Each set was 
assigned to a particular site based on the known anisotropy of the quadrupole splitting 
[in] . In figure El Hi is plotted against the angle 9 between /fext and the c-axis at 
several temperatures. By symmetry, H2{0)=Hi{—9) and H^{6)=H4{6) (not shown). 
The ^-dependence can be well fitted by a simple sinusoidal curve (curves in figure 12)), 

Hi + H2 = Hu cos 26 + const , 

Hi~ H2 = Hssm2e . (4) 

T-dependences of Hu and Hg are plotted in figure El 

The 6'-dependence of Hi is not symmetric around ^ = 0, leading to a finite Hg. Both 
mf'^ and m^'^ contribute to H^ through the uniform (diagonal) and the staggered 
(off-diagonal) parts of the coupling Aij, respectively. If m^'^ and mf'^ had the same 
T-dependence, this should stand for H^ and H^. This is indeed not the case, since Hg 
extrapolates to a large finite value at T=0, while Hu almost vanishes similarly to the 
magnetic susceptibility [4 . This proves that a sizable staggered moment mf'^ is induced 
at T=0 at field values H^xt <^ H^ while the uniform moment is nearly zero. This can be 
qualitatively explained by considering Hi as a perturbation for an isolated dimer [T7j. 
Since Tii has matrix elements between singlet and triplets on the same dimer, triplets 
are mixed into the ground state, resulting in finite mf'^ to first order in perturbation. 
In contrast, m^'^ is much smaller than mf'^, since only higher order terms contribute 
to m^'^. Although field-induced staggered moments have been reported in the Haldane 
chain compound NENP IT^ , this is the first observation in quasi 2D systems. 

For a quantitative analysis, we write 

Hi-H2 = 2 {B^, + C^,) {mt sin 9 + m^^ cos 6) 

+ 2 {B^^ - a^) sin e + 2 {B,, - C,,) m^^ cos 9 

+ AD^z (m^^ sin 9 + mf^ cos 9) + AD-jcyVn^y sin 9, (5) 

where we have defined B=Aii, C=Ai2, and D=Ai3 and used symmetry properties to 
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Figure 3. Temperature dependences of Hg, Q iJsu, and i/^s (jSJ at -ffext=6.948 T. 



obtain other Aj^'s. The sin26' dependence of H1-H2 in Q is reproduced if the moments 
vary as m^, m^, m;^^'^ oc cos6', and m^^, m;^^^ oc sin 6'. This is approximately confirmed 
by the numerical calculations presented below. We define 

a, ^ mt{e = 0) = ml{e = 0) , a, = -mi{e = n/2) , 
M. ^ mtf{e = 0) , ^ mti^\d = n/2) . 

Combining these relations with (0), we obtain 

Hg = Hsu + Hss (6) 
Hsu = (S,. + a. + 2I^..) + (S,, + a,) M,^ + 2D,,Mf (7) 

-f^ss = {Bxx ~ Cxx + '^Dxy) CTx ~^ {Bzz — Czz) O'z- (8) 

The analysis of the ^^B NMR spectrum in the 1/8-plateau phase [20] indicates that 
all hyperfine couplings except B are approximately given by the classical dipolar fields 
that we can calculate. Since the diagonal components of B+C+2D were determined 
previously all the coupling parameters in the above equations are known except 
for Bxz- At high temperatures, the staggered moments ax, (Jz are much smaller than 
the uniform moments M^, M^'^, and . Assuming ax=az=Q above T=70 K, 

which is also supported by preliminary numerical calculations, we determined Bxz from 
the data of Hs and the magnetization at 70 K, using ((7j). The parameters in ((7j) and (jHl) 
are determined as Bxz+Cxz+'2.Dxz=^-'^^, S^.2+C^2=0.09, 2L'^^=0.01, Bxx-Cxx+'2Dxy= - 
0.06, and i?2^-C'2^=-0.32 (T/yUs)- We then obtained Hsu at other temperatures (figureEI) 
from the magnetization data. Finally, subtracting this from Hs-, Hss is determined in 
the whole temperature range as shown in figure El This demonstrates that the staggered 
magnetization increases very steeply below 10 K and saturates at lower temperatures. 

Measurements of Hi were extended to higher fields up to 26 T. The magnetization 
data in ref. [S] was used to calculate Hsu |2I|; which was then subtracted from the data 
of Hs to obtain Hss- In the main panel of figure 01 we plot Hsu and Hss in the low 
temperature limit. 

In order to see if these results are explained by anisotropic interactions, we have 
calculated rrij by exact diagonalization (ED) of the hamiltonian Hq+Hi+H' at T=0 




Figure 4. Field dependences of Hg^ and Hss measured at low temperatures, 2.5 K 
(for iJoxt=7 and 12 T), 1.5 K (15 and 17.9 T), 0.64 K (20 T), 0.46 K (24 T), and 
0.21 K (26 T). The curves are obtained from the ED results of a^, cr^, M^, and M^^^ 
shown in the inset by using and (©. 



and -ffext=6.948 T in the (llO)-plane for clusters with up to 24 sites for various values 
of D. Values of other parameters are fixed: J=85 K, J'=54 K 0, L'7J=-0.02 [12], and 
5's=0.023. By requiring that the calculated value of Hg from (0) to (jSI) agrees with the 
experimental data at the lowest temperature (2.5 K), we obtained -D/J=0.034. Taking 
into account the uncertainty in Qg leads to possible values of D/J between 0.030 (for 
gs = 0.030) and 0.038 (for gs = 0.014). The choice of D' has little effect on the calculated 
values of rrij. 

Using the same parameter values, we have also calculated rrii at higher fields, as 
shown in the inset of in figure ^ The agreement between the calculated results of Hsu 
and Hgs shown by the lines in the main panel of figure |3] and the experimental data 
is satisfactory up to about 20 T. Thus we now have a quantitative explanation for the 
development of a uniform magnetization at fields as low as 18 T. Note that the staggered 
moment is sizable even at low fields (~ 0.03 /is/Cu at 15 T, figure lU inset), where the 
uniform moment is negligible. The experimental and theoretical results of Hgs deviate 
at higher fields, probably due to finite size effects as discussed below. 

Let us now discuss the dynamics. The nuclear spin-lattice relaxation rate (1/Ti) was 
measured for i/ext || c up to 27.5 T. It shows an activated behavior at low temperatures 
(figure El inset). The activation gap (A) shown in figure El follows a linear if-dependence 
below 20 T, A = Aq — g^fisH with Ao=34.8 K and gz=2.28 as expected. At higher 
fields, however, it deviates from this relation and remains finite even above the expected 
critical field Hc=Ao/ gzfJ^B='22.7 T. The gap reaches a minimum at 26.5 T, which is at 
the boundary to the 1/8-plateau, and increases again inside the 1/8-plateau (27.5 T). 
The finite gap above He, which naturally explains the absence of Bose condensation, is 




Figure 5. Field dependence of the excitation gap (solid circles in the main panel) 
determined from the activated behavior of l/Ti (inset) are compared with the results 
of ED calculation (open circles). The data at 27.5 T, where NMR spectrum show many 
lines in the 1/8-plateau phase, were obtained for the line with the largest (negative) 
hyperfine field. 



consistent with earlier reports. For example, the gap was estimated from the specific 
heat data [HI as A=4.6 (3.2) K at i/ext=22 (24) T. The ESR data HO] also shows a 
similar deviation from the linear if-dependence of the lowest triplet energy near He. 

One expects that the singlet-triplet mixing due to Tii will prevent the gap from 
closing, as has been discussed for NENP We have confirmed this by numerical 

calculations. The energy of the lowest excited state obtained by ED calculations using 
the same parameters is plotted in figure |31 The constant value at low fields is due to the 
singlet bound state of two triplets having lower energy than one triplet 23^. However, 
they do not contribute to the nuclear relaxation. Near 22 T, the results clearly indicate 
level repulsion due to mixing, reproducing the experimental behavior. Let us note one 
subtle point. The interdimer DM interaction Ti' splits the one-triplet excitations into 
two branches Only one of them is mixed with the singlet by Tii and this branch 
must have a lower energy for the level repulsion to occur. This determines the sign of 
D' to be negative. Note that the sign of D' was not known before. 

Finally, the theory is not applicable when the field is too close to the 1/8-plateau. 
The cluster size of our calculation is not large enough to take proper account of 
interaction between triplets, which should become important near the 1/8-plateau. 

In conclusion, we have shown that a large staggered magnetization is induced by 
the magnetic field in the presence of the intradimer DM interaction and staggered g- 
tensor, whose magnitudes are obtained as 0.03 < D/J < 0.038 for 0.014 < Qs < 0.03. 
The singlet-triplet mixing and level repulsion caused by these interactions account for 
the finite uniform magnetization below He and the persistence of the excitation gap 
above H^. They thus turn the quantum phase transition at He into a crossover. This 
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is the first example of sucli field-induced phenomena in quasi 2D spin systems. Our 
determination of the magnitudes and sign of the DM interactions should set the stage 
for further investigation of other remarkable properties of SrCu2(B03)2. In particular, 
the signature of a phase transition in the specific heat observed above the 1/8-plateau 
by Tsujii et al. [TTj calls for a precise understanding of the interplay between the 
DM interaction and other bosonic interaction such as bound state formation and 
correlated hopping Plj. 
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